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Abstract 

Working within the framework of the QCD light cone sum rules (LCSR), we 
compute and discuss the nonfactorizable higher twist effect in B° — > D*°tt° to 
make an all-around examination of its role in the charmed color-suppressed B 
decays. Analogously to the case of B° —* D°ir , such effect turns out to be of 
the same strong phase as the factorizable amplitude, and modifies constructively 
the magnitude by order (40 — 90)% so that the effective coefficient a{ = C\ + 
C2/3 receives a positive correction comparable numerically with it. Nonleading 
as the soft effect in question is, our findings for it, along with the previous LCSR 
analyses of B° — > D°ir°, are suggestive of the dominance of soft exchanges in these 
charmed color suppressed B decays. Also, the emphases are put on importance of 
understanding intensively various related higher twist and transverse momentum 
effects to interpret the data onB^ D 0( * )(7r°, 77, 77'), J/ip K^*h 
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The formulation of the QCD factorization (QCDF) pQ and soft-collinear effective 
theory (SCET) |2j has greatly renewed and deepened our understanding of QCD dy- 
namics in the two-body hadronic decays of B-mesons, and as such considerably enhanced 
our confidence in precisely estimating the decay amplitudes and CP-asymmetry effects. 
Especially, ones have essentially changed minds in understanding of the charmed color- 
suppressed B decays and realized that the nonfactorizable soft exchange contributions 
dominate over the hard-exchange ones. For all these progresses, however, being unable 
to deal with non-perturbative QCD dynamics from the first principle makes yet it dif- 
ficult to theoretically interpret the experimental observations of the color suppressed 
process B° -> D 0( -*°h° One attempts to make sense of them by a global analysis of 
the data on B — > D^ix decays, either in a model-independent manner [3J E] , or in 
the framework of SCET [7j and Perturbative QCD [S]. The leading soft effects in such 
decays are supposed to include final-state rescattering processes and nonfactorizable con- 
tributions induced by the color-octet operators. Resorting to various phenomenological 
approaches, such as one-particle exchange model |5] and Regge theory [TU], one models 
the soft contributions from the former. QCD sum rule approach is initially used ^Tj and 
QCD light-cone sum rule (LRSR) method [12] is subsequently applied [13] to obtain an 
estimate of the latter. Focusing on the B° — > D°tt° case, the author of ^Tj finds that 
the soft exchanges between the color-octet quark pair cq and the pion, which concern 
higher-twist components of the pion, provide a negative correction at a level of 70% 
of the factorizable amplitude, while the LCSR approach predicts a positive numerical 
effect at (50 — 110)% order JT3]. As compared to the B — > J/^/K case where the effect 
in question contributes to the magnitude by (30 — 70)% [T^j, the higher- twist effect in 
B°^D° 7T° is, as it were, more important, as expected. Also, the similar calculations are 
carried out for B — > titt [T3*l E3 and B — > ^7] , but the resulting numerical impacts 
are found to be far less than those in both B° — > D q tt° and B — > J/^/K result 
in accordance with the argument for QCD factorization. 

Whereas the nonfactorizable higher-twist effect, as estimated, plays a role comparable 
numerically with the factorizable amplitudes of the color-suppressed B° — > D°7r° and 
B — > J/^/K decays, in the letter we would like to examine how it affect B° — > D*°7r°, to 
get an all-around understanding of its role in the charmed color suppressed B decays 

We start with looking briefly back at the dynamical features of the color suppressed 
B° — > D°(*°^7r . In the Fierz transformed form, the relevant effective weak Hamiltonian 
reads [IB] 
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- — f= V cb V* d 



(Cx(/0 + \c 2 {h))om + 2CM6M] , (i) 
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where are the Wilson coefficients, the CKM matrix elements, and 

o 1 = (cr^)(Jr^), 6, = (c^r^)(d^b), (2) 

the color-single and -octet quark operators, respectively, with T M = 7 M (1 — 75) and A a 
being the color SU(3) matrices. It is clear that the transitions induced by 0\ come from 
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gluon exchanges and thus the corresponding matrix elements are non-factorizable. Since 
the ejected _D°(°*) in the decays is a heavy meson inside which the two quarks are asym- 
metric in momentum distribution, the non-factorizable vertex diagrams prove themselves 
to be QCD infrared divergent so that QCD factorization breaks down 19 . This hints 
at the dominance of soft dynamics. If we introduce an effective coefficient a 2 such that 
the decay amplitude of B° — > D*°ir is made be of the following parametrization: 

A(B°(p + q) - D*°(e,p)7r°(q)) = -GpV^e ■ qm D ,f D * F*^{m D *)a 2 , (3) 

with e, m£>* and fp* being the polarization vector, mass and decay constant of D* meson, 
respectively, and F B ^ n (m 2 D *) the B — > 71 form factor defined as 

2 2 

nrri fry* « 

(Tr(g) I u^b I B(p+q)) = (2q+p),F 1 B ^(p 2 )+p, B ^ " (F^tf) - Ff^tf)), (4) 

a 2 is actually supposed to absorb all the nonfactorizable contributions including higher 
twist component we are to calculate. 

Let's go over to LCSR calculation of the nonfactorizable higher-twist effect involved 
in B° — > D*°tt°. With this end in view, we opt for the following correlator: 

U^{p,q,k) = i 2 J J d i xd' l ye- l{p+q)x+l(p - k)y 

x (7r°(q)\T{u(y)^c(y), O x (0), m b b(x)t l5 d(x)}\0) , (5) 

where the current operators «7 M c and bi^^d are used to interpolate D*° and B° meson 
fields respectively, and an unphysical momentum k is introduced to make B channel free 
of pollution by D*ir resonances, but it must disappear from the physical B° — > D*°7r° 
amplitude in the dispersion integral. Schematically, the higher twist contributions to the 
correlator can be illustrated by Feynman diagrams in Fig.l. The correlation function 
I1 M is a functions of 6 independent invariants, which, adapting to the present purpose, 
are to be chosen as p 2 , q 2 , (p + q) 2 , (p — k) 2 , k 2 and P 2 = (p + q — k) 2 . At small light- 
cone distances x 2 ~ y 2 ~ (x — y) 2 ~ 0, or equivalently in large and space-like regions 
of P 2 , (p + q) 2 and (p — k) 2 , the Operator Product Expansion (OPE) is applicable to 
Eq.(|5|). Furthermore, not violating the validity of OPE we simply let k 2 be zero, and 
take chiral limit approximation q 2 = to the pion. 

Inserting on the right hand side of © a complete set of hadronic states with D*° 
quantum numbers, we obtain the phenomenological expression U^ D *\ 

7, P^Ks,(p + q) 2 ,P 2 ,p 2 ) 

H(D*) ^ ' 

S 
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where 

Gf, = i^e x I d 4 xe-^ x (D*°(e x ,p-k) 7r (g)|T{6i(0), m b b(x)i l5 d(x)}\0) 
x J 

= n(( P + ?) 2 ,P 2 ,p 2 ,e*w-g < ) (7) 

A 

is a two point function with Q { = p, q and k, p^ £> *- ) and s^ 15 ^ are respectively 
the hadronic spectral function and threshold mass squared of the higher states in D* 
channel. The integral of (|B|L invoking quark-hadron duality, is supposed to coincide with 
the s > Sq part of the QCD spectral representation of I1 M 

1 7 lm s U^ CD \s,(p + q) 2 ,P 2 ,p 2 ) 

n r s-ip-w ' (8) 



S D* 



with s®* being the QCD effective threshold in D* channel, and ^lm s Il^ CD ^ the related 
QCD spectral density. Then we equate © with the QCD answer for (j5J) and make the 
Borel transformation (p — k) 2 — > T, yielding 

G,(p + q) 2 ,P 2 ,p 2 ) = - JL ^- T: ds e^-^ T Irn s Ilf CD \s,(p + q) 2 ,P 2 ,p 2 ). (9) 

To proceed, saturating the two-point function G M with the intermediate states of 
B° quantum numbers, we have its hadronic expression in which the matrix element 
(D*°(e x ,p — k)n (q)\Oi\B (p + q)) enters the lowest pole term. When the condition 
(p+q) 2 = (p + q — k) 2 = m 2 B is imposed on the B°(p+q) — ► D*°(e x ,p — k)n°(q) transition, 
the unphysical momentum k vanishes automatically so that the corresponding matrix el- 
ement is made physical. On some algebraic manipulations we see easily that the required 
physical amplitude can be parameterized in term of a function factor M^~ k ' in front of 
(p — k)^ in the kinematic decomposition of e ^(D* (e x ,p — k)ir (q)\Oi\B°(p + q)). This 

denotes that we need only to calculate the related part ^qcd °f n^ CD \ Furthermore, 
to apply quark hadron duality to B channel we stick down Im. s Il^ CD \s, (p + q) 2 , P 2 ,p 2 ) 
in a form of dispersion integral in variable (p + q) 2 , and make the analytic continuation 
P 2 — > m B , keeping (p + q) 2 fixed and letting p 2 = m 2 D *. Then using the standard 
procedure the desired sum rule can be achieved, if the relevant QCD double spectral 
density Pqcd{s, s',m 2 B ,p 2 ) is known, where s' is a spectral variable in B channel. 

Now we embark on the computation of Hqcd t° extract Pqcd- It is advisable that we 
work with the fixed-point gauge in which the light cone expansion of quark propagator 
|2Uj . which takes in a correction term due to the interaction with one background-field 
gluon. After some lengthy calculations we obtain the twist-3 contribution: 

n^> = -^h- r ds ^ r\ dt r 1 du r ^ 



,w:i 8y/2n* J ml s - (p - k) 2 J 4 J x (s,t,ml) J x (s,t, m %) [m 2 - {p + qu) 2 }v 
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mz. 



X 03,(1 - u, u - v, v) I (s - - 1) + 2t[(p + qY - v 



x (2v-x(s,t,m%))} 



(10) 



In the above equation, f 3n is a nonperturbative quantity denned by the vacuum-pion 
matrix element of the local operator ucr^j^Gapd, 3vr refers to a pionic twist-3 dis- 
tribution amplitude related to the nonlocal operator d{Qi)a^^G a p{vy)u{x) |2T], and 

2 

X = { s ~ ^r)/{ s ~ m s)- The fact that the threshold Sq is far below the squared mass 
of b quark allows us to make an expansion of (10) in \. Then passing over those terms 
vanishing after Borel transformations we have, up to order 0(x 3 ), 



p tw3 (s,s',m B ,p ) 



1 f u dv 
x<2/ —(p 3 „(l-u,u-v,v) 
Jo tr 



m b f; 



16v^7T 2 (s' -J9 2 ) -I- 



dt 



mz. 



2vt(s' - p 2 



+2 



u dv 



(1 - 3t)(s' -p ) / —4> 3n (l-u } u-v } v) 



o v 



mz 



+ (t - 1) [S ^03tt(1 -u,u-v,v] 



v=0 



x(s,t,m 2 B ) 



+ 



m 2 \ d 



(t - 1) i s - -f i — (^jMi 



+2(2t - l)(s' - p 2 )^0 3 .(l - u, w - v, v) 



X' 



(11) 



v=0 



with it = (m 2 . —p 2 )/{s' — p 2 ). Following the same line the twist-4 contribution is attain- 
able. Here we do not give it any more to save some space. 

At the end, the required nonfactorizable effect in B° — > D* 7t° can be numerically 
estimated by the following LCSR answer: 



(D* (e,p)7r (q)\d 1 (0)\B (p + q )) 



e ■ q mi, mn* 



1 ^ e (m|-K-m^»(l-n))/n)/T' 

M ( f U 



4:V2TT 2 f D *f B m 2 B (m 2 B - m 2 D ,) 
dtl f 3n 



ds e (.m» D .-8)/T 



— 4> 3n (l -u,u-v,v) 
v z 



mz — m 



x 



D 



-[2vt+(l-3t) X ]+[s-^\ (1-t) 



+(t-l) 



mz 



s — 



^3tt 



(1 — u,u — v, v) 



v=0 
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( 771^ 777^ X 

+ 2 (2t - 1)-* — ^<Ml - «, « - v, v) 



u w- 



1 

— 3 tt(1 - U, W- V,V) 
V 



X 



v=0 



2 2(2t 



u dv ( 
o i> V 



(2t — l)0j_(l — W, M — U, f ) + 0J_(1 — 1i, U — f , f ) 



(1 — u, u — v, v) + 0_|_(1 — m, u — v, v) x 



1 / 3 r 



+ 2 V 

l_d 

v dv 



(2t — l)(f)±(l — u, u — v , v ) + — u, u — v, v) 
(2t — l)<f>±(l — u, u — v, v) — 0J_(1 — u, u — V, V 



X 



v=0 



(12) 



where Uq = (m 2 — m? D ,)/(sQ — m 2 D S) with Sq being the QCD effective threshold in B 
channel, T' indicates the Borel parameter corresponding to the variable (p + q) 2 , the 
pionic distribution amplitude <p± (<j>±) is of twist-4 [2*T] . 

To go on with numerical discussions. We take [221 123 m D = 2.007 GeV, m c = 
1.3 ± 0.1 GeV, f D . = 240 ± 10 MeV , s° = 6 ± 1 GeV 2 and T = 1 - 2 GeV 2 for the 
D*° channel parameters, and m B = 5.28 GeV, m b = 4.7 ±0.1 GeV, f B = 180 ±30 GeV, 
s B = 35 ± 2 GeV 2 and X" = 8 — 12 GeV 2 for the B channel ones. The nonperturbative 
quantities concerning % meson are chosen as |2H]: fn = 132 MeV, fz-ni^b) = 0.0026 GeV , 
and 



V?37r(ai,/ib) = 360aia 2 a;3 



l + ^(7«3-3) 



<P±(ati, n h ) = 305 2 («i - a 2 )a 2 [- + 2e 7r (l - 2a 3 )], 



<fx(ai,fM b ) = 30<5 2 (1 - a 3 )a 2 



26* (1 - 2a 3 ) 



(13) 



whth 



??? 



B 



hit, 



m b /2, a scale characteristic of the typical virtuality of the 



underlying b quark, and cu 37T = —0.28, 5 2 = 0.17GeV 2 , e w = 0.36. The other parameters 
are fixed as C^) = -0.26, C 2 (fx b ) = 1.12 and F^{m 2 D ,) = 0.30 [221 - 

Using these inputs, we see that with both the Borel parameters varying in the re- 
spective intervals, Eq. (jl2j) reveals a good stability in numerical value. From the sum 
rule windows, the resulting impact on the decay amplitude could be derived using 



As{ B° - /TV) = V2G F V cb V: d C 2 (D*°(e,p)n (q)\d 1 (0)\B°(p + q)). 



(14) 



Nevertheless, for the purpose of comparing with the factorizable amplitude and reducing 
the uncertainty, it is more advisable to calculate the ratio of As over the factorizable 



5 



amplitude Af- The result is 

R = A S (B° -> D*\°)/A F (B° -> ,0*%°). 

= 0.43 - 0.95. (15) 

The range quoted here is due to the variations of the various inputs and is obtained by 
adding linearly up the corresponding uncertainties. The soft contribution is shown to be 
of the same phase as the factorizable one, indicating that it does not provide the decay 
amplitude with any strong phase, and with a certain uncertainty such effect modifies the 
magnitude of the factorizable amplitude at a comparable or even identical level. 

When a comparison is drawn with the situation of B° — > D°n° in which the ratio 
is estimated at between 0.54 and 1.15 [TJ], the decay in question receives a little less 
soft modification within the framework of LCSR, as expected. Also, It is easily figured 
out by intuitive picture that a bit more ratio is observed in the above two cases than 
R = 0.30 - 0.70 yielded for B -> J/ipK in [TS]. 

To more give prominence to the role of the soft effects, we might, as in Eq(3), express 
A = A F + A s as 

A(B° -> /TV) = ~G F V cb V* d (ml - m^/i^^K^, (16) 

in term of a phenomenological parameter = a{ + 2C 2 R 2 , with a{ = C x + C 2 /3 
corresponding to the factorizable contribution and R 2 being introduced to parameterize 
the soft effect. While R 2 turns out numerical smallness, the resulting influence on a 2 is 
considerable, due to the large coefficient 2C 2 multiplying it. Numerically, the soft effect 
is such that it adds a positive number to a 2 , changing the value of a 2 from 0.12 to 
0.16 — 0.23. Such a result, along with that obtained in B° — > D°tt , keeps still away from 
that demanded by the experiment data | a 2 (A(B° -> D°(*°h°) \= 0.57 ± 0.06. 

At this point, a minor comments should be in order on a couple of open problems 
with soft dynamics in charmed color suppressed B decays: As a result of the dominance 
of nonfactorizable soft dynamics in B° — > D°(*°^ir , in addition to transverse momentum 
dependence which is important to get a result of infrared finiteness, the higher-twist 
components have to be counted in the light cone wavefunction description of _D°(*°) 
meson. The same must also be done for B — > J/ipK, notwithstanding the fact that no 
corresponding infra divergence appears so as to validate QCD factorization for it, owing 
to the emitted quarkonium system (CC) being of a small size ~ 1/ 'ot s {m c )m c -C 1/Aqcd- 
On the other hand, as shown in B — ► -kit, J/ipK, the hard spectator contributions have 
a divergent integral from the end point region of the distribution amplitude for the light 
mesons, when the twist-3 components is included. However, the transverse momenta 
may not be omitted in the region. We believe that for B° -> D ^ 7r it is this case too. 
The other resources of soft effects, of course, such as final state interactions [24J, have to 
enter consideration in an effective manner for understanding the data. 

Colligating the LCSR analyses made for B° — > D ^* ^ , J/ip K, we conclude that 
nonfactorizable higher twist corrections, in spite of being a nonleading effect, are im- 
portant to understand the charmed color suppressed B decays. These findings favor, 
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in particular, the argument for QCD factorization: lange-distance QCD dynamics pre- 
dominates in the color-suppressed processes such as B° — > D 0( -*°\tt , rj, rj'). A further 
investigation is necessary to improve the present results. The staple improvement is to 
take into account the QCD radiative correction to the correlator. The other mends are 
going to wait until the higher-twist amplitudes of pion become updated and available. 
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Figure 1: Feynman diagram illustrating nonfactorizable higher twist contributions to 
the correlator (6). Solid lines represent quarks, dashed line gluon, wavy lines stand for 
external currents. The square refers to the insertion of the operator Oi, and oval the 
pionic higher twist wavef unctions. The cross denotes the other point to which gluon line 
can be attached. 
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